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(Dated: February 7, 2008) 

Correlations in the azimuthal angle between the two largest transverse momentum jets have been 
measured using the D0 detector in pp collisions at a center-of-mass energy •Js = 1.96 TeV. The 
analysis is based on an inclusive dijet event sample in the central rapidity region corresponding 
to an integrated luminosity of 150 pb~^. Azimuthal correlations are stronger at larger transverse 
momenta. These are well-described in perturbative QCD at next-to-leading order in the strong 
coupling constant, except at large azimuthal differences where soft effects are significant. 
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Multi-parton radiation is one of the more com- 
plex aspects of perturbative Quantum Chromodynamics 
(pQCD) and its theory and phenomenology are being ac- 
tively studied for the physics programs at the Fermilab 
Tevatron Collider and the CERN LHC 01 • The proper 
description of radiative processes is crucial for a wide 
range of precision measurements as well as for searches 
for new physical phenomena where the influence of QCD 
radiation is unavoidable. A clean and simple way to 
study radiative processes is to examine their impact on 
angular distributions. Dijet production in hadron-hadron 
collisions, in the absence of radiative effects, results in 
two jets with equal transverse momenta with respect 
to the beam axis (pt) and correlated azimuthal angles 
dijet = |<^jeti — 0jet2| = TT. Additional soft radiation 
causes small azimuthal decorrelations, whereas A0 dijet 
significantly lower than tt is evidence of additional hard 
radiation with high pT- Exclusive three-jet production 
populates 27r/3 < Aipdijct < ti" while smaller values of 
dijet require additional radiation such as a fourth jet 
in an event. Distributions in A(/)dijGt provide an ideal 
testing ground for higher-order pQCD predictions with- 
out requiring the reconstruction of additional jets and 
offer a way to examine the transition between soft and 
hard QCD processes based on a single observable. 

A new measurement of azimuthal decorrelations be- 
tween jets produced at high pT in pp collisions is pre- 
sented in this Letter. Jets are defined using an iter- 
ative seed-based cone algorithm (including mid-points) 
with radius 7?.conc — 0.7 0. The same jet algorithm 
is used for partons in the pQCD calculations, final-state 
particles in the Monte Carlo event generators, and recon- 
structed energy depositions in the experiment. A(/) dijet 
is reconstructed from the two jets with highest pr in 
an event. The observable is defined as the differential 
dijet cross section in A(/) dijet, normalized by the dijet 
cross section integrated over A(/) dijet in the same phase 
space (l/(7dijet) (rfcrdijet/rfA(^ dijet). (Theoretical and ex- 
perimental uncertainties are reduced in this construc- 
tion.) Calculations of three-jet observables at next-to- 
leading order (NLO) in the strong coupling constant ctg, 
have recently become available Quantitative com- 

parisons with data yield information on the validity of 
the pQCD description and increase sensitivity for gaug- 
ing potential departures that could signal the presence of 
new physical phenomena. 

Data were obtained with the D0 detector ^ in Run II 
of the Fermilab Tevatron Collider using pp collisions at 
^/s = 1.96 TeV. The primary tool for jet detection was a 
compensating, finely segmented, liquid-argon and ura- 
nium calorimeter that provided nearly full solid-angle 
coverage. Calorimeter cells were grouped into projective 
towers focused on the nominal interaction point for trig- 
ger and reconstruction purposes. Events were acquired 
using multiple-stage inclusive-jet triggers. Four analysis 
regions were defined based on the jet with largest pr in 




dijet (rad) 

FIG. 1: The dijet distributions in four regions of p™'"^. 
Data and predictions with p™^'' > 100 GeV are scaled by suc- 
cessive factors of 20 for purposes of presentation. The solid 
(dashed) lines show the NLO (LO) pQCD predictions. 



an event {p^^^) with the requirement that the trigger 
efficiency be at least 99%. The accumulated integrated 
luminosities for events with p^'^^ > 75, 100, 130, and 
180 GeV were 1.1, 21, 90, and ISOpb^^ (±6.5%), re- 
spectively. The second leading px jet in each event was 
required to have pt > 40 GeV and both jets were re- 
quired to have central rapidities with |?/jet| < 0.5 where 
yjct = 5 In {{E + p^)/ [E - p^)) and E and p^ are the 
energy and the longitudinal momentum of the jet. 

The position of the pp interaction was reconstructed 
using a tracking system consisting of silicon microstrip 
detectors and scintillating fibers located within a 2 T 
solenoidal magnet. The vertex z-position was required to 
be within 50 cm of the detector center which preserved 
the projective nature of the calorimeter towers. The sys- 
tematic uncertainty associated with the vertex selection 
efficiency is less than 3% for At/) dijet > 27r/3 and « 8% 
for A(/) dijet ~ 7r/2. The missing transverse energy was 
calculated from the vector sum of the individual trans- 
verse energies in calorimeter cells. Background from cos- 
mic rays and incorrectly vertexed events was eliminated 
by requiring this missing transverse energy to be below 
Q.l p^^^. Background introduced by electrons, photons, 
and detector noise that mimicked jets was eliminated 
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FIG. 2: Ratios of data to the NLO pQCD calculation for 
different regions of p™^^- Theoretical uncertainties due to 
variation of fj.r and /i/ are shown as the shaded regions; the 
uncertainty due to the PDFs is indicated by the solid lines. 
The points at large dijct are excluded because the calcula- 
tion has non-physical behavior near the divergence at tt. 



based on characteristics of shower development expected 
for genuine jets. The overall selection efficiency is typi- 
cally ~ 83% for A(/)dijot < 57r/6 and drops to « 78% as 

A^dijot 71". 

The pt of each jet was corrected for calorimeter show- 
ering effects, overlaps due to multiple interactions and 
event pile-up, calorimeter noise effects, and the energy 
response of the calorimeter. The calorimeter response 
was measured from the pt imbalance in photon -I- jet 
events. The relative uncertainty on the jet energy cali- 
bration is w 7% for jets with 20 < pr < 250 GeV. The 
sensitivity of the measurement to this calibration was re- 
duced by normalizing the dijct distribution to the inte- 
grated dijet cross section. Nevertheless, this provides the 
largest contribution to the systematic uncertainty (< 7% 
for A(/) dijct > 47r/5 but up to 23% for Ac/) dijct < 27r/3). 

The correction for migrations between bins due to fi- 
nite energy and position resolution was determined from 
events generated with the herwig and PYTHIA 
programs. The generated jets were smeared according 
to detector resolutions @. The angular jet resolution 
was determined from a full simulation of the D0 detec- 
tor response, ft was found to be better than 20 mrad 



for jets with energies above 80 GeV. The jet pr res- 
olution was measured from the pt imbalance in dijet 
events. It decreases from 18% at px = 40 GeV to 9% for 
Pt = 200 GeV. Finite jet pt resolution can lead to am- 
biguities in the selection of the two leading pT jets. This 
effect is large at small A0 dijct where contributions from 
higher jet multiplicities dominate. The generated events 
were reweighted to describe the observed A0 dijct distri- 
bution. This provided a good description of the observed 
Pt spectra of the four leading pT jets. The correction for 
migrations is typically less than 8% for Ai/i dijet > 27r/3 
and w 40% for A0 dijct ~ 7r/2 with a model depen- 
dence of less than 2%. Only for pf'''' < 130 GeV and 
at A0 dijet ~ 7i'/2, is the model dependence as large as 
» 14%. 

The corrected data are presented in Fig. Q as a func- 
tion of A^ dijct in four ranges of p™'''^. The inner error 
bars represent the statistical uncertainties and the outer 
error bars correspond to the quadratic sum of the sta- 
tistical and systematic uncertainties. The spectra are 
strongly peaked at A(/) dijet ~ tt; the peaks are narrower 
at larger values of p'^^^. Overlaid on the data points 
in Fig. n] are the results of pQCD calculations obtained 
using the parton-level event generator nlojet-) — h 
and CTEQ6.1M ^ parton distribution functions (PDFs) 
with as{Mz) = 0.118. The observable was calculated 
from the ratio of the predictions for 2 — + 3 processes 
(rfo-dijet/rfA0 dijct) and 2-^2 processes (cTdijct), both at 
leading order (LO) or NLO, 
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CTdijct 



rfcTdijet 



(N)LO 



dA(/) dijet 



(N)LO 



The renormalization and factorization scales are chosen 
to he fir = IJ-f = 0-5 p^'^^. The ratio is insensitive to 
hadronization corrections and the underlying event [loj . 

NLO pQCD provides a good description of the data. 
As shown in Fig. [21 data and NLO agree within 5- 
10%. The theoretical uncertainty due to the PDFs 
is estimated to be below 20%. Also shown is the ef- 
fect of renormalization and factorization scale variation 
{0.25 p^^^ < fir J < Pt'^^)- The large scale dependence 
for A(/) dijet < 27r/3 occurs because the NLO calculation 
only receives contributions from tree-level four-parton fi- 
nal states in this regime. Results from pQCD at large 
A0 dijet in Figs, n and [5] were excluded because fixed- 
order perturbation theory fails to describe the data in 
the region A(j)dijet ~ tt where soft processes dominate. 

Monte Carlo event generators, such as herwig and 
PYTHIA, use 2^2 LO pQCD matrix elements with phe- 
nomenological parton-shower models to simulate higher 
order QCD effects. Results from herwig (version 6.505) 
and PYTHIA (version 6.225), both using default param- 
eters and the CTEQ6L PDFs, are compared to the 
data in Fig. |3| herwig describes the data well over the 
entire A0 dijet range including A0 dijct ~ t^- pythia with 
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FIG. 3: The A(^dijct distributions in different p™^" ranges. 
Results from herwig and pythia are overlaid on the data. 
Data and predictions with p™^" > 100 GeV are scaled by suc- 
cessive factors of 20 for purposes of presentation. 

default parameters describes the data poorly — the distri- 
bution is too narrowly peaked at A^^ijct ~ and lies sig- 
nificantly below the data over most of the A^dijot range. 
The maximum pT in the initial-state parton shower is 
directly related to the maximum virtuality that can be 
adjusted in pythia. The shaded bands in Fig.|31indicate 
the range of variation when the maximum allowed virtu- 
ality is smoothly increased from the current default by a 
factor of four llj. These variations result in significant 
changes in the low A^dijot region clearly demonstrating 
the sensitivity of this measurement. Consequently, global 
efforts to tune Monte Carlo event generators should ben- 
efit from including our data. 

To summarize, we have measured the dijet azimuthal 
decorrelation in different ranges of leading jet pt and 
observe an increased decorrelation towards smaller pr- 
NLO pQCD describes the data except for very large 
Ac/) dijet where the calculation is not predictive. 
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